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ABSTRACT 



We investigate magnetization as functions of external magnetic field H in 
the [/-infinite Anderson lattice model within the leading order approximation 
in the 1/iV-expansion. At T = 0, at 7J = Hm where the Zeeman energy is 
equal to a certain characteristic energy in the system, the magnetization curve 
has a kink and the differential susceptibility dM/AH shows a jump. At finite 
temperature, dM/d-ff shows a peak around Hyi- Its maximum value increases 
with decreasing T and saturates to a finite value at T ^ 0. When H < Hm, 
the / and the conduction electrons form the rcnormalized bands with a large 
Fermi surface determined by the Luttinger sum rule. On the other hand, when 
H > Hm, the bands reform themselves significantly free from the Luttinger 
sum rule, eventually leading to a small Fermi surface aX H Hm- The results 
are consistent with the metamagnetic properties observed in the heavy fermion 
CeRuaSis. 
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There have been extensive studies on the metamagnetic behavior observed in 
the normal metalhc states of the heavy fermion compounds. The investigations 
have been mostly focused on CeRu2Si2,||l|, ||, ||, Q whose magnetization shows 
an abrupt increase with increasing iJ at a certain critical value, H = Hm ~ 
7.7T. Many other properties of the material also show pronounced changes at 
H ~ Hm- The drastic is the change of the Fermi surface observed in the recent 
de Haas- van Alphen (dHvA) experiments: Q the large Fermi surface observed at 
H < Hm seems to suggest existence of the itinerant /-electrons, while the small 
Fermi surface observed sd, H > Hm seems to suggest the localized /-electrons. 
Several theoretical studies have been already done on these problems. [|[ ||, 0, ^ 
^, [l^ They comment on the magnetization curves but nothing on the drastic 
change of the Fermi surface. 

The purpose of the present paper is to investigate theoretically the meta- 
magnetic properties including both the magnetization and the Fermi surface. 
We employ the U - infinite Anderson lattice model in the auxiliary boson 
representation. ||l^ Previously, the model has been studied in the limit of zero 
external magnetic field within the approximation including terms of the leading 
order in the expansion with respect to the inverse of the spin-orbit degeneracy 
(1/A'^-expansion) under the strict local constraints guaranteeing the equivalence 
of the bosonic version to the original [/-infinite model. The results very well 
account for the various properties of the heavy fermion systems in the high tem- 
perature incoherent regime exhibiting the Kondo effects, the low temperature 
coherent regime with heavy fermions and the intermediate temperature regime 
showing the crossover behavior between the two limiting regimes in a unified 
way. The present study is a straightforward extension of the previous study to 
cases with finite external magnetic field. 

Our Hamiltonian is given by 

,7 J 

m——J km 'i m— — J 

J 

+^l'E E E(^^™^"''"''''<j-^^+'^-^-), (1) 

i m——J km 

where m stands for the degrees of freedom due to the spin-orbital degeneracy, 
m = — J, — J+1,--- J (2 J -1-1 = N), and i stands for the lattice sites, 
i = 1,2,- •• A^L- fcm represents the generalized wave vector specifying wave 
vector as well as spin-orbital degrees of freedom, is the creation operator 
for the c-electron. bf and f^-^ are the creation operators for the slave-boson (SB) 
and the pseudo-fermion (PF) representing the empty and the singly occupied 
states of the i-th /-site, respectively. Here we set g/^B = 1 and ffcMB — 0. Then, 
mH is the Zeeman energy for the /-electron under the external magnetic field 
H.^ The energies e^^ and e / are measured relative to the chemical potential 

/X. 
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This model (j^) is equivalent to the original {/-infinite Anderson lattice model 
as long as it is treated within the physical space where the following local con- 
straints hold: 

J 

Qr= + ^i^^ = = 1' 2, • • . TVl). (2) 

m— — J 

The expectation value of an operator O calculated under the local constraint 
(I) is given by0 

(0)= hm {d\lQ.,)x/{\{Q.)x, (3) 

i I 

where {A}x = Tr[e-'3^M]/Tr[e~^^^] with Hx = H + E., ^^Q^■ For calculat- 
ing {A)\, we treat 14^ as a perturbation and employ the standard perturba- 
tion method using the Feynman diagrams together with the expansion {1/N- 
expansion) from the large limit of the spin-orbit degeneracy N, while keeping the 
total degrees of freedom for the c-electrons to be constant: Sfc„, 1 = 2. 

Following the procedure mentioned above we calculate the single particle Green's 
functions within the leading order in power of The explicit forms of the 

c-electron, the SB and the PF Green's functions are given, respectively, byp^ 

Gk^{iujn) = [iujn - ek,„ -T,k^{iujn)]~^, (4) 
B{ivn-Xt) = [iM„ - Aj - ^n™(ij^„ - Aj)]"\ (5) 

m 

Frn{iuJn ~ Xi) = [lUJn - - + TOiJ]~\ (6) 

where the self-energy parts are given by 

X ^^;„(iu;„ -f ii^„' - \^)B{iiyn, - \i)/{Q^)x\, (7) 

with {Qi)x = {btbi)x + Y.Tn{finJim)>.- The occupation numbers of the SB and 
the PF are given in the leading order of the 1 /A^-expansion by 

{bth)x = -r^S(ii.„-A,), (9) 

+ r^i?(i.„-A.) ^""'j;"";^-^ (10) 

a(ll'n) 
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Note that Fm is not modified in the lowest order of l/N as shown in eq.(^, 
while the correction of 0{1/N) to {fj^fim)x must be included in eq.(p^) to 
obtain {Qi)\ of 0((1/-/V)°). Equations fl)-([lO|) constitute a set of self-consistent 
equations (SCE), which was already solved in the limit oi H = analytically at 
T = and numerically at finite temperatures. Now, we solve the SCE with 

To begin with, we investigate the magnetization process at T = 0. With 
-ff 7^ 0, we rewrite the occupation numbers, eqs.(^ and (p^), and the sum of 
them, {Qi)\, in the low temperature limit as 

{fLhn)x c-'='(^'+^/-"^)+An/„e-'^(^'+^/-^°), 

m 

(11) 

with Aufm = —a ■^Kellm{uj)\^_^^_^^, where Eq and a are, respectively, the 
binding energy and the residue of the resonance state in the SB spectrum, 
ImB(i/ ~ Xi + iO+), which are determined by the relations: 

ef-Eo-Y] ReU^ief - Eo) = 0, (12) 



V:^Ren™(c.) 



(13) 



In the above, both Eq and a are functions of T and H to be determined later. 
Eq = Eq\t=h=o represents the characteristic energy of the system, which cor- 
responds to the Kondo temperature defined in the single impurity Anderson 
model. We note that the continuum in the SB spectrum, which has finite 
intensity for v > Xi + ef — min{EQ, JH), is irrelevant to calculate the occupation 
numbers at T = 0. From eq.(^), the average number of the /-electron with m-th 
orbit is given by njm = liTai),.^ao[{f^fim)x/{Qi)x]- Using eq.(|ll|), we obtain 
at T = as 

Arif 
1 



/m = -TT ^n^fjSr,i,J, (14) 



where n^j = a/(l -f a) is the incoherent part of the /-electron with m = J, 
which is due to the zeroth order term w.r.t. V^^, and a = limr^o c'^''^^^^"-'. 
Equation ( [l4| ) yields the average number, nj = n fm , and the magnetization, 
M — mnf,n, of the /-electron, respectively, at T = 0: 

Uf = il-a)/{l + a)+n%, (15) 

M = ^mAn/,„/(l + a) + Jn'^j. (16) 
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\Vk,\^B{ef -JH- iLJn)njj6.mj, (17) 



Substituting eqs.(||), (^) and (|T^) into eq.(|^), we find 

at T = 0, with Eq^ = Eq — mH. Substituting eq.(p7|) to eq.(||), we obtain the 
c-electron Green's function. 

Now we solve cqs.(p^) and ( p^ ) with eq.(^ at T = to obtain Ea and a 
together with a for given H. Hereafter for the numerical calculations we assume 
simple fcm-dependence as: e^^ = Cc — and |Vfc^p — V^^k„, together with a 
square density of states: p°„(0 = J2k^ ^i^ - ^/c„) = for < C < 1, 
and Pm(^) = for otherwise. The parameters are chosen as: J — ^ {N = 
6),n = n^ + Hf = 1.5, ec - £/ = 0.9 and V = 0.2.0 

In Fig.l, we show the self-consistent solutions at T = for a and a as 
functions of H, and those for Eo by the dashed line in Fig. 2. We find that, at 
H < Hm, Eq > JH and a = 0, while, at > Hm, Eq = JH and a ^ 0, 
where Hu = 11.12 x lO"'^ is of order of E°/2. {E° = 22.77 x lO"'' as plotted 
in Fig. 2.) Each of a, a and Eq has a kink a,t H = Hm- Correspondingly, the 
magnetization curve has a kink at iJ = Hm as shown in Fig. 3 (see the dashed 
line for T — 0). At H = Hm, the slope of M increases abruptly with increasing 
H, because so does that of the incoherent part n'jj (or a) as shown in Fig.l 
contributing to the 2nd term in eq.(|l6|). 

Next, we investigate the magnetization process at finite temperature. In this 
case, we can not use the analytic expressions eqs.(^^-(|l7|) obtained at T = 0. 
Then we solve the SCE (|)-(|l^ numerically in the same way as taken in the 
previous studies. Fig. 2 shows iJ-dependence of Eq calculated from eq.(^2|) 
at several temperatures. When H > Hm, we observe that {Eq — JH) oc T at the 
low temperatures, and find that the extrapolated value of ePiJ^~^o) at T — > 
coincides with a calculated through the SCE at T = (see also Fig.l). In Fig. 3, 
the magnetization curves are plotted at several temperatures. At the lower 
temperatures, the slope of M increases significantly around H = Hm, while, 
at the higher temperatures, M increases only monotonically. In Fig. 4, we plot 
the differential susceptibility defined by dM/dH as functions of H at several 
temperatures. At the lower temperature T < 5 x lO^'', dA//diJ shows a peak 
around H = Hm- Its maximum value increases with decreasing T and saturates 
to a finite value at T ^ 0. Such behaviour in dM/dH at T — > has been 
observed in the recent magnetization measurements at very low temperature. Q 

The magnetic susceptibility Xs is plotted as function of T in Fig. 5, where 
Xs = J{J + l)n'j/3T is the contribution due to the incoherent part of the /- 
electron, rij — "/t?x' while Axs is the coherent part due to the lowest order 
corrections in the 1/A'^-expansion. We observe a maximum in Xs at T = Tmax) 
where T^ax ~ 5.5 x lO"'' is of order of E^/A. At T < r^ax, ^Xs is the 
dominant term leading to the enhanced Pauli paramagnetism, Axs ^ l/-^0' 
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while, at T ;2> Tmax, Xs is the dominant term leading to the Curie law, Xs ~ 
J(J + l)/3r.|Q At T ^ Tmax, the incoherent part, n°, contributes dominantly 
to M even for H ^ 0. Thus, the metamagnetic behaviour due to the rapid 
increase of n'jj is no more distinguished at T > Tmax, as seen in Fig. 4 and 
as observed in the experiments. Q We emphasize that both Hm and T,„ax are 
proportional to Eq which determine the Fermi liquid properties of the system, 
e.g., Xs|t=o ~ l/-E'o- Such relationship among Hm, T^i^^ and the characteristic 
energy of the Fermi liquid has been observed in the pressure dependence of 
metamagnetic properties.]^ 

Finally, we discuss the quasi-particle properties and the Fermi surface. Sub- 
stituting eq. (|l^) into eq. (||) , we obtain the c-electron Green's function at T = 
for H < Hm as Gk,„ {iuJn) = L-y^i ^1 P'' 



E2„ 



with 



E 



± _ 



Eom)/{E^ 



EL 



In Fig. 6(a), E^ are plotted as functions of for H < Hm, which describe 
the coherently c/-hybridized excitation with the Zeeman split resonance level 
Eom- On the other hand, for H > Hm, the imaginary part of the c-electron 
self-energy with m = J is finite at w < as shown in eq.(pT|) (note that the 
continuum in the SB spectra is finite for > Xi + Cf —min{Eo, JH)). Therefore 
the c/-hybridized excitation with m = J becomes incoherent. However the 
incoherent excitation has a finite energy gap as seen in Fig. 6(b), because the 
resonance level Eqj — 0. Thus they are irrelevant to the low energy properties 
at T = 0. 

Using eqs.(|^) and (|l^), the average number of the /-electron with m-th orbit 
is written by: 



T'fm 



dSfc^(iu;„) 
d(icj„) 



Then the total number of the c- and the /-electrons per site is given by: 



T 



dlo g{Gfc„^(ic^„)-i} 
d(it^„) 



When T ^ and H < Hm, it yields + nf ^ N'^ J^m Efc„ ^(-^kj^ by 
using iTy — and ImEfc^^(iO+) ~ 0. Thus the Fermi surface is determined by 
the total number of the c- and the /-electrons, leading to a large Fermi surface 
as seen in Fig. 6(a). This is nothing but the Luttinger sum rule.[^ On the 
other hand, when T = and H > Hm, n° = n°j =/= 0. Therefore, the Luttinger 
sum rule breaks down. As nP^j increases with increasing H, the phase volume 
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enclosed by the Fermi surface with J decreases abruptly, leading to a small 
Fermi surface determined almost only by the c-electrons as seen in Fig. 6(b). The 
drastic change of the Fermi surface around H = Hm has been observed in the 
recent dHvA experiments.]^ 

In conclusion, we have investigated the magnetization process of the Ander- 
son lattice model in the metallic phase within the leading order approximation 
in the 1/A^-expansion. A critical field Hm has been shown to exist at T = 0: 
a.t H > Hm, the incoherent part of the /-electron due to the zeroth order term 
w.r.t. the c/-mixing becomes relevant to the physical quantities, although it 
is irrelevant at H < Hm- Consequently, the magnetization increases far more 
steeply with increasing H iov H > Hm, and the differential susceptibility shows 
a jump at iJ = Hm- Moreover the phase volume enclosed by the Fermi surface 
has been found to change drastically at H — Hm- at H < Hm, it is as large as 
determined by the Luttinger sum rule, while, at H > Hm, it shrinks leading to 
a small Fermi surface determined almost only by the c-electrons. At finite tem- 
perature, the smooth metamagnetic behavior has been observed at T < Tmax, 
the temperature where the magnetic susceptibility has a maximum. The results 
agree with the observation in the metamagnetic heavy fermion CeRu2Si2. 
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Figure Captions 

Fig.l Self-consistent solutions for the residue a of the resonance level of the 
slave-boson (dotted line) and a defined in eq.(|l4|) (solid line) at T = as 
functions of H. 

Fig. 2 Self-consistent solutions for the binding energy Eo of the resonance level 
of the slave-boson at T = (dotted line) and at several finite temperatures 
(solid lines) as functions of H. 

Fig. 3 Magnetization M at T = (dotted line) and at several finite tempera- 
tures (solid lines) as functions of H. 

Fig. 4 Differential susceptibility dM/dH at T = (dotted line) and at several 
finite temperatures (solid lines) as functions of H. 

Fig. 5 Temperature dependence of the magnetic susceptibility. The dashed line 
is the contribution of the incoherent part of the /-electron, x^, the dotted 
line is that of the coherent part of 0((l/iV)°), Axs, and the solid line is 
the sum of them, Xs- 

Fig. 6 Schematic structures of the renormalized bands E^ for m = —J,—J + 
1, • ■ • J at H < Hm (a) and H > Hm (b) as functions of S^k- The dashed 
lines in (b) are for m — J. The bare c-electron band is also plotted 
(dotted line). 
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